Fast domain wall propagation under an optimal field pulse in magnetic nanowires 
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We investigate field-driven domain wall (DW) propagation in magnetic nanowires in the framework 
of the Landau-Lifshitz- Gilbert equation. We propose a new strategy to speed up the DW motion 
in a uniaxial magnetic nanowire by using an optimal space-dependent field pulse synchronized with 
the DW propagation. Depending on the damping parameter, the DW velocity can be increased by 
about two orders of magnitude compared to the standard case of a static uniform field. Moreover, 
under the optimal field pulse, the change in total magnetic energy in the nanowire is proportional 
to the DW velocity, implying that rapid energy release is essential for fast DW propagation. 

PACS numbers: 75.60.Jk, 75.75.-c, 85.70.Ay 



Recently the study of domain wall (DW) motion in 
magnetic nanowires has attracted a great deal of atten- 
tion, inspired both by fundamental interest in nanomag- 
nctism as well as potential industrial applications. Many 
interesting applications like memory bitsQ, or mag- 
netic logic devices [3[ involve fast manipulation of DW 
structures, i.e. a large magnetization reversal speed. 

In general, the motion of a DW can be driven by a mag- 
netic field 0-@| and/or a spin-polarized current 0- 11 1. 



Although the DW dynamics in systems of higher spatial 
dimension can be very complicated, some simple but im- 
portant results were obtained by Schryer and Walker for 
effectively one-dimensional (ID) situations [l^]: At low 
field (or current density), the DW velocity v is linear in 
the field strength H until H reaches a so-called Walker 
breakdown field H w 12]. Within this linear regime, DW 
propagates as a rigid object. For H > H w , the DW 
loses its rigidity and develops a complex time-dependent 
internal structure. The velocity can even oscillate with 
time due to the "breathing" of the DW width. The time- 
averaged velocity v decreases with the increase of H, re- 
sulting in a negative differential mobility, v can be again 
linear with H approximately when H 3> H w . The pre- 
dicted v-H characteristic is in a good agreement with 
experimental results on permalloy nanowires [4- 6]. Re- 
cently a general definition of the DW velocity proper for 
any types of DW dynamics has been also introduced [l3j]. 

For a single-domain magnetic nanoparticle (called 
Stoner particle) , an appropriate time-dependent but spa- 
tially homogeneous field pulse can substantially lower 
the switching field and increase the reversal speed 
since it acts as an energy source enabling to overcome 
the energy barrier for switching the spatially constant 
magnetization [hH 15 1. In the present letter, we inves- 
tigate the dynamics of a DW in a magnetic nanowire 
under a field pulse depending both on time and space. 
As a result, such a pulse, synchronized with the DW 
propagation, can dramatically increase the DW velocity 
by typically two orders compared with the situation of a 
constant field. Moreover, the total magnetic energy typi- 
cally decreases with a rate being proportional to the DW 



velocity, i.e. the external field source can even absorb 
energy from the nanowire. 
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FIG. 1: A schematic diagram of two dynamically equivalent 
ID magnetic nanowire structures. (A) Easy axis is along the 
wire axis (z-axis); (B) Easy axis (z-axis) _L the wire axis (x- 
axis). The region between two dashed lines denotes the DW 
region. 

A magnetic nanowire can be described as an effectively 
ID continuum of magnetic moments along the wire axis 
direction. Magnetic domains are formed due to the com- 
petition between the anisotropic magnetic energy and the 
exchange interaction among adjacent magnetic moments. 
Let us first concentrate on the case of a uniaxial mag- 
netic anisotropy: Two dynamically equivalent configura- 
tions of ID uniaxial magnetic nanowires are schemati- 
cally shown in Fig. [TJ Type A shows the wire axis to 
be also the easy-axis (z-axis). Type B shows the easy 
axis (z-axis) is perpendicular to the wire axis (x-axis). 
Although our results described below apply to both con- 
figurations, we will focus in the following on type B. 
The spatio-temporal dynamics of the magnetization den- 
sity M(x,t) is gove rned by the Landau-Lifshitz-Gilbert 
(LLG) equation [ 
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where [7] = 2.21 x 10 5 (rad/s)/ (A/ m) the gyromag- 
netic ratio, a the Gilbert damping coefficient, and M s 
is the saturation magnetization density. The total ef- 
fective field H t is given by the variational derivative of 
the total energy with respect to magnetization, H t = 



2 



— (5E/5M)/ fig, where no the vacuum permeability. The 
total energy E = J_ dxe(x) can be written as an inte- 
gral over an energy density (per unit section- area), 
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where x is the spatial variable in the wire direc- 
tion. Here K, J are the coefficients of energetic 
anisotropy and exchange interaction, respectively, and 
H is the external magnetic field. Moreover, we have 
adopted the usual spherical coordinates, M(x,t) = 
M s (sin 9 cos (f>, sin 9 sin 0, cos 9) where the polar angle 
9(x, t) and the azimuthal angle 4>(x, t) depend on position 
and time. 

Hence, the total field H t consists of three parts: the 
external field H, the intrinsic uniaxial field along the 
easy axis H K = (2KM Z / ijlq)z, and the exchange field 
H J which reads in spherical coordinates as 12|, ll7| , 
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Following Ref. [12], let us focus on DW structures ful- 
filling Ocf) /Ox — 0, i.e. all the magnetic moments rotate 
around the easy axis synchronously. Then the dynamical 
equations take the form 



T9 = a H, 
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where we have defined T = (1 + a 2 )|7| _1 , and Hi(i = 
r, 9, <p) are the three components of the external field 
in spherical coordinates. In the absence of an exter- 
nal field, an exact solution for a static DW is given by 
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exp(cc/A) where A = ^Jj/(KM 2 ) is the width 
of the DW. We note that a static DW can exist in a con- 
stant field only if the field component along the easy axis 
is zero, H z = 0. In fact, according to Eqs. (|4j) static solu- 
tions need to fulfill = [implying <p — ta,n~ 1 (H y /H x ) 
is spatially constant] and 
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or, upon integration, 
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cos 29 + H r (9) = constant. (6) 



Considering the two boundaries at 9 = 0(x — > — oo) and 
9 = tt(x — > +oo) for the DW, we conclude H r (0) = 



H r (ir), which requires H z = 0. In this case, the station- 
ary DW solutions under a transverse field are described 

asi = J[^J{KM 2 sin 2 9 - fi Q M s H sin9)/ J}- 1 d9. 

Thus, when an external field with a component along 
the easy axis is applied to the nanowire, the DW is ex- 
pected to move. We use a travelling-wave ansatz to de- 
scribe rigid DW motion [l2^. 
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where the DW velocity v is assumed to be constant. Sub- 
stituting this trial function into Eq. (j4]), the dynamic 
equations become 

TsmOv = -A(aH e T sin 9<j) — aH^ — H g . (8) 

Eq. © describes the dependence of the linear velocity v 
and the angular velocity tf> on the external field H. Our 
following results discussion will be based on Eqs. ((5J. 

Let us first turn to the case of a static field case applied 
along the easy axis (z-axis in type B of Fig. [TJ, Hg = 
—Hsm9, = 0. Here we recover the well-known static 
solution for a uniaxial anisotropy [3], 
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where the azimuthal angle <f>(t) — (f)(0) + \j\Ht/ (1 + a 2 ) is 
spatially constant (i.e. 0cf>/0x — 0) and increases linearly 
with time. 

Let us now allow the applied external field to depend 
both on space and time. Our task is to design, under a 
fixed field magnitude H, an optimal field configuration 
H(x,t) to increase the DW velocity as much as possible. 
From Eqs. ijHJ), we find a manifold of solutions of specific 
space-time field configurations described by a parameter 



H r (x,t) =Hcos9 1 H 8 (x,t) = -H sin0/Vl + u 2 , 

= -uiTsin0/Vl + u 2 . (10) 

The velocities v and <f> reads 
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The previous static field case is recovered for u = 0. The 
maximum of the velocity v m with regard to u us reached 
for u = 1/a, 
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where the angular velocity is zero, cA = 0. On the other 
hand, <j) attains a maximum for u = —a, where, in turn, 
the linear velocity vanishes. In Fig. [5] we have plotted 
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FIG. 2: (Color online) The DW propagation velocity v ver- 
sus the parameter u at the different damping values a — 
0.1,0.2,0.5,0.8. The other parameters are chosen as A = 
20nm and H = lOOOe. 



the dependence of the velocity on the parameter u for 
different damping strengths and typical values for the 
DW width A and the magnitude H of the external field. 

To understand the physical meaning of the maximum 
velocity v m , we note that, according to Eqs. ©, the held 
components Hg and are required to be proportional 
to sin# to ensure the constant velocity under the rigid 
DW approximation. Moreover, at u = 1/a we have Hg = 
aHtf,, and from the identity 

(aH e + H^) 2 + (aH^-Hg) 2 = (1 + a 2 )(H 2 - H 2 ), (13) 

we conclude that the term (aHg + H^) is maximal re- 
sulting in a maximal velocity according to Eqs. ((5J. As 
a result, the new velocity under the optimal field pulse is 
larger by a factor of v m /v = y/l + a 2 / a « 1/a compared 
to a constant field with the same field magnitude. To 
give a practical example, the typical value for the damp- 
ing parameter in permalloy is a = 0.01 which results in 
an increase of the DW velocity by a factor of 100. 

It is instructive to also analyze the optimal field pulse 
according to Eq. (fTQ|) with u = 1/a in its cartesian com- 
ponents, 



H x (x,t) = iTsin20(l - a/y/1 + a 2 )/2, 
H y (x,t) = -Hsin6/y/l + a 2 , 



(14) 



H z (x,t) = H (cos 2 9 + a sin 2 0/y/l + a 2 ), 

where 9 follows the wave- like motion tan ^fe^l = e xp(^ — 
Ji+a 2 ^' ^ n El we plotted these quantities at t = 
around the DW center where the main spatial variation 
of the pulse occurs. Note that the space-dependent field 
distribution should move with the same speed v m syn- 
chronized with the DW propagation. Near the DW center 
the components H x and H z are (almost) zero whereas a 
large transverse component H y is required to achieve fast 
DW propagation. Qualitatively speaking, the transverse 
field causes a precession of the magnetization resulting in 



its reversal. This finding is consistent with recent micro- 
magnetic simulations showing that the DW velocity can 
be largely increased by applying an additional transverse 
fieldGl. 
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FIG. 3: (Color online) The x,y,z components of the optimal 
field pulse. The parameters are chosen as a — 0.01, A = 
20nm for permalloy [f|. The field magnitude is H = lOOOe. 

It is also interesting to study the energy variation un- 
der the optimal field pulse, 



dE 
~dt 




H t +M-—\=P a +P h . 



(15) 

The first term P a is the intrinsic damping power due to 
all kinds of damping mechanisms described by the phe- 
nomenological parameter a. According to the LLG equa- 
tion P a = J+^d^dM^ ig alwayg negativeO, 

implying an energy loss. Ph is the external power due to 
the time-dependent external field. From Eq. (fTTj). both 
powers are obtained as 
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such that the total energy change rate is 
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Note that the intrinsic damping power is independent of 
the parameter u and always negative, whereas the total 
energy change rate is proportional to the negative DW 
velocity. Thus, for positive velocities (u > —a) the total 
magnetic energy decreases while it grows for negative ve- 
locities (u < —a). In the former case energy is absorbed 
by the external field source while in the latter case the 
field source provides energy to the system. The optimal 
field source helps to rapidly release or gain magnetic en- 
ergy which is essential for fast DW motion. This aspect is 
very different from the reversal of a Stoner particle where 
the time-dependent field is always needed to provide en- 
ergy to the system to overcome the energy barrier [l4j|. 
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Moreover, our new strategy of employing space- 
dependent field pulses can also be applied to uniaxial 
anisotropies of arbitrary type: Let w(9) be the uniaxial 
magnetic energy density. The static DW solution in the 
absence of an external field reads x = J x 1 (9)d9, where 



X (0) = y/[w(6)-w ]/J. 



(19) 



Here wo is the minimum energy density for magneti- 
zation along the easy axis. By performing analogous 
steps as before, we obtain the the optimal velocity as 

v rn = J 7 !/* ! where Xmax denotes the maximum of 

X(0) throughout all 9. 

On the other hand, our approach is not straightfor- 
wardly extended to the case of a magnetic wire with biax- 
ial anisotropy. To see this, consider, a biaxial anisotropy 
£i = -KM? + K'M? where the coefficients K, K' cor- 



respond to the easy and hard axis, respectively [12|. The 
LLG equations read 
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Let us assume 4>(x,t) 
mined by the applied field. Substituting the travelling- 
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, where now A 
20l we obtain 



Tsin9v = -A(aHg + H, 
olK'M s sin 6* sin 20o /mo ■+ 
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(21) 

(aff# - Hg) = 0. (22) 



For a static field along z-axis Hg = — H sin 9, H<j, = 0, the 
solution is just the Walker's result v — \j\AH/a (Note 
here A also depends on -ff)[l2|. To implement our new 
strategy, we need to find the maximum of the right-hand 
side of Eq. (f2Tj) under two constraints of Eq. (|22|) and 
Eq. (fT3|) with Hg and H$ being proportional to sin 6*. 
The unique solution to this problem is indeed a constant 
field along the z-axis which is thus the optimal field con- 
figuration. 

In summary, our theory is general and can be applied 
to a magnetic nanowire with a uniaxial anisotropy which 
can be from shape, magneto-crystalline or the dipolar 
interaction. The experimental challenge of our proposal 
is obviously the generation of a field pulse focused on 
the DW region and synchronized with its motion. How- 
ever, the field source synchronization velocity can be pre- 
calculated from the material parameters. As for the re- 
quired localized field (See Fig. [3]), we propose to employ 



a ferromagnetic scanning tunneling microscope (STM) 
tip to produce a localized field perpendicular to the wire 
axis[20| and use a localized current to produce an Oer- 
sted field along the wire axis 21|. Moreover, such re- 
quired localized fields may also be produced by nano- 
ferromagnets with strong ferromagnetic (or antiferromag- 
netic) coupling to the nanowire. We also point out that, 
although the field source typically does not consume en- 
ergy but gain energy from the magnetic nanowire, the 
pulse source may still require excess energy to overcome 
effects such as defects pinning, which is not included in 
our model. At last, the generalization of the strategy be- 
yond the rigid DW approximation, and to DW motion 
induced by spin-polarized current will also be attractive 
direction of future research. 
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